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ABSTRACT 

Wc present a numerical study of the properties of the stellar velocity distribution in 
stellar discs which have developed a saturated, two-armed spiral structure. We follow 
the growth of the spiral structure deeply into the non-linear regime by solving the 
Boltzmann moment equations up to second order. By adopting the thin-disc approx- 
imation, we restrict our study of the stellar velocity distribution to the plane of the 
stellar disc. We find that the outer (convex) edges of stellar spiral arms are character- 
ized by peculiar properties of the stellar velocity ellipsoids, which make them distinct 
from most other galactic regions. In particular, the ratio ci : <J2 of the smallest versus 
largest principal axes of the stellar velocity ellipsoid can become abnormally small (as 
compared to the rest of the disc) near the outer edges of spiral arms. Moreover, the 
epicycle approximation fails to reproduce the ratio a^^ : arr of the tangential versus 
radial velocity dispersions in these regions. These peculiar properties of the stellar 
velocity distribution are caused by large-scale non-circular motions of stars, which in 
turn are triggered by the non-axisymmetric gravitational field of stellar spiral arms. 

The magnitude of the vertex deviation appears to correlate globally with the 
amplitude of the spiral stellar density perturbations. However, locally there is no simple 
correlation between the vertex deviation and the density perturbations. In particular, 
the local vertex deviation does not correlate with the local gravitational potential and 
shows only a weak correlation with the local gravitational potential gradient. We find 
that the local vertex deviation correlates best with the spatial gradients of mean stellar 
velocities, in particular with the radial gradient of the mean radial velocity. 

Key words: Galaxies: general - galaxies: evolution - galaxies: spiral - galaxies: 
kinematics and dynamics 



1 INTRODUCTION 

Already since the end of the nineteenth century studies of 
stellar kinematics in the solar neighbourhood revealed that 
the velocity dist ributio n of stars in galactic discs is non- 
isotropic fKobold ll890l ). An ea rly interpretation by Kapteyn 
(|l905h and Eddington (| 19061 ') invoked a superposition of 
(isotropic) stellar streams with different mean velocities, by 
this creating an anisotropic velocity distribution. How ever, 
an alternative interpretation by Schwarzschild (| 19071 ) be- 
came the general framework for describing the local velocity 
distribution of stars of equal age, though a contamination 
by stellar streams is still considered to be important for a 
proper analysis of velocity data. The Schwarzschild distribu- 
tion is based on a single but anisotropic ellipsoidal distribu- 



tion function. It is characterized by a Gaussian distribution 
in all three directions U (radial r), V (tangential (j)) and 
W (vertical z) in velocity space, but with different veloc- 
ity dispersions arr, a<t,<t>, and azz- In general, the velocity 
distribution is described by a velocity dispersion tensor 

O-lj = (Vi - Vi){vj - Vj), (1) 

where i and j denote the different coordinate directions, i.e. 
radial (r), azimuthal (0) and vertical (2) direction and Vi 
gives the corresponding velocity. The bar denotes local av- 
eraging over velocity spac^^e emphasize that our definition 
of the velocity dispersion tensor is differe nt from the defi- 
nition of e.g. Kuijken & Tremaine (|l994 ) who do not use 
a square mark. Hence, our a^j is equivalent to Kuijken & 
Tremaine's aij. . The principal axes of this tensor form an 
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imaginary surface that is called the velocity ellipsoid. The 
principal axes of the velocity ellipsoid need not to be aligned 
with the coordinate axes. In that case, the vertex deviation 
Iv (which is defined as the angle between the direction in 
velocity space pointing from the Sun to the Galactic centre 
and the direction of the aij major principal axis) does not 
vanish. 

The general trend for the velocity dispersions is that 
the radial velocity dispersions are the largest and the verti- 
cal velocity dispersions are the smallest among stars of same 
age and stellar type. This trend holds when dispersions are 
averaged locally , i.e. in the solar neighbourhood (Dehnen 
& Binnev [l99a) . or globally over the entire galactic discs. 
For instance, the ratio arr '■ cr^^ : a^z of the radial, tan- 
gential, and vertical stellar velocity dispersions in the disc 
of N GC 488 is approximately 1 : 0.8 : 0.7 (Gerssen et al. 
1 19971 ). 

In the case of stationary, axisymmetric systems and 
appropriate distribution function^ (DF) the velocity ellip- 
soid is aligned with the coordinate axes. If the stellar or- 
bits are nearly circular (i.e. with small radial amplitudes 
with respect to the guiding centre), then the Oort ratio 

= (700 : a^^ of tangential versus radial stellar velocity 
dispersions can be calculated in zero order approximation 
within the epicyclic theory to be X^p — —B/{A — B), where 
A and B are the usual Oort constants. For a flat rotation 
curve one gets Xfp = 1/2. More realistic axisymmetric dis- 
tribution functions give first order corrections, which yield 
values of X' ^ = 0. 59 (or 0.66) for the solar vicinity (Kuijken 
& Tremaine ll99ll . hereafter KT91). Clearly, these values ex- 
ceed those predicted by the epicycle approximation. On the 
other hand, observations in the Milky Way (MW) give lower 
(mean) values around X'^ — 0.42 — 0.45, which are substan- 
tiall y belo w the epicycle value = 1 /2 (Kerr & Lynden- 
Bell Il986l : Ratnatunga & Upgren Il997l ). The exact reason 
for this discrepancy is still unknown. It might stem from a 
DF asymmetric in (KT91) or from invalid assumptions, 
i.e. a violation of stationarity or axisymmetry. 

Additionally, non-vanishing vertex deviations were 
found i n the solar vicinity in many stu dies (e .g. Wielen [l974l . 
Mayor Il972l . Ratnatunga & Upgren 1 19971 '). These vertex 
deviations could be explained by spiral structures (Mayor 
Il970l . Yuan Il97lh . by this supporting the importance of 
non-axisymmetry of the underlying gravitational potential. 
Moreover, a clear relation between 1^ and arr has been 
found. For smaller velocity dispersions larger vertex devi- 
ations have been measured, e.g. for a radial velocity dis- 
persion of about 15 km s~^ one gets a vertex deviation of 
25° -30° (Wielenll974|). 

A major step in analysing the local stellar velocit y spac e 
was the astrometric satellite mission Hipparcos (ESA ll997h . 
Earlier results concerning the general behaviour of the ve- 
locity ellipsoid, e.g. the dependence of the dispersions or 
the vertex deviations as a func tion of B — V w ere corrobo- 
rated (Dehnen & Binney Il998l . Bienayme Il999l . Hogg et al. 
I2OO5I ). Moreover, the more numerous known proper motions 



^ i.e. those DFs only depending on the two integrals of motion 
energy E and the 2-component Lz of the angular momentum or 
those obeying a third integral of motion which is symmetric under 
a simultaneous change of sign of Vr and Vz 



allowed for a detailed mapping of the velocity space in the 
solar vicinity. It turned out that the local velocity distribu- 
tion of stars exhi bits a lot of substructure (Dehnen Il998l . 
Alcobe & Cubarsi lioosi ). For example, two major peaks in 
the velocity distribution were found where the smaller sec- 
ondary peak is well detached by at least 30 kms ~^ from the 
main peak (the "u anomaly"). Dehnen (|l999l ) attributed 
this bimodality to the perturbation exerted by the Galactic 
bar assuming that its outer Lindblad resonance (OLR) is 
close to the Sun. Miihlbauer & Dehnen l|2003f ) showed that 
a central bar can also explain vertex deviations of about 10° , 
and Oort ratios below 1/2. Recent investigations about the 
influence of the Galactic bar on Oort's C constant corrob- 
orated these results and gave stronger limits on the bar's 
properties (Minchev et al. 120071 '). Minchev & Quillen (|2007t ) 
showed also that a stationary spiral cannot reproduce the "u 
anomaly". Though less likely, an alternative interpretation 
of the bimodality as a result of non-stationary spiral arms 
cannot be ruled out completely. 

In general, a non-axisymmetric gravitational potential 
might lead to both, a misalignment of the velocity ellipsoid 
and an additional correction of the velocity dispersion ratio 
X^ with respect to the value predicted by the epicycle ap- 
proximation. In the case of spi ral pe rturbations, this conclu- 
sion was made e.g. by Mayor l| 19701 . for a review of different 
mechanisms see KT91) and numerically confirmed recently 
by Vorobyov & Theis (|2006l . hereafter VT06). Using linear 
perturbation analysis KT91 studied different discs subject 
to oscillating perturbations (like a rotating bar or a spiral) 
or subject to a static non-axisymmetric dark matter halo. In 
the case of large-scale oscillating perturbations the predicted 
vertex deviations oscillate, too. Once the axisymmetric back- 
ground gravitational potential and the pattern speed of the 
perturbation are fixed, its amplitude depends in a complex 
way only on the perturbed mean velocity and its spatial 
gradient. For a static elliptic potential Kuijken & Tremaine 
(l994h showed that the observed mean kinematic data {X 
and Zv) can be explained, if the ellipticity of the disc is about 
q ~ 0.88 and the Sun's location is near to the minor axis of 
the corresponding potential. It is remarkable that the correc- 
tion terms introduced by the non-axisymmetry are negative 
(i.e. in agreement with the observations in the Milky Way), 
whereas the higher order terms derived from a more realis- 
tic ax isymmetric DF are positive. Similarly, Blitz & Spergel 
(|l99ll ) argued for a slowly rotating mildly triaxial Galactic 
potential, resulting in an outward motion of 14 kms~^ of 
the LSR and a vertex deviation of 10°. 

However, it is not clear if the non-axisymmetric grav- 
itational field is entirely responsible for the observed ver- 
tex deviation. The existence of moving groups of stars was 
shown to p roduce large vertex deviations (Binney & Mer- 
rifield Il998l ) . The situation may become even more compli- 
cated because moving groups of stars may in turn be caused 
by the non-axisymmetric gravitational field of spiral arms. 
Therefore, a detailed numerical study of the vertex deviation 
in spiral galaxies is necessary. 

In this paper, we present the stellar hydrodynamics sim- 
ulations of linear and non-linear stages of the spiral structure 
formation in stellar discs by solving numerically the Boltz- 
mann moment equations up to second order in the thin-disc 
approximation. In contrast to linear stability analyses, we 
follow the growth of the spiral structure self-consistently. 
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Compared to A''-body simulations, we have the advantage 
that we can follow the spiral structure formation from a very 
low perturbation level deeply into the non-linear regime. The 
former is usually hardly accessible by self- consistent A''-body 
simulations due to particle noise. We study the local evolu- 
tion of the velocity dispersion tensor in the plane of the disc, 
which includes the anisotropy (X^), the ratio of the smallest 
versus largest principal axes of the stellar velocity ellipsoid, 
and the vertex deviation (Zv). 

The paper is organized as follows. In Section [2] we give 
a short description of the numerical code. The model galaxy 
and initial conditions are described in Section |3l The phys- 
ical explanation for the growth of the spiral structure in 
our model disc is given in Section |4l The shape and orien- 
tation of stellar velocity ellipsoids is studied in Section [5l 
In Section [6] we discuss the applicability of the epicycle ap- 
proximation and possible causes of the vertex deviation in 
spiral galaxies. The summary and conclusions are given in 
Section [3 



2 THE METHOD AND NUMERICAL CODE 

Stars are coUisionless objects that move on orbits deter- 
mined by the large-scale gravitational potential. An exact 
description of such a system requires the solution of the col- 
lisionless Boltzmann equation for the distribution function 
of stars f{r,v,t) in the phase space {r,v). A general solu- 
tion of the Boltzmann equation, which is defined in a six- 
dimensional position and momentum phase space, is usually 
prohibited due to the large computational expense. How- 
ever, taking moments of the Boltzmann equation in velocity 
space turns out to yield a set of equations that is numerically 
tractable. 

In our recent paper (VT06) we presented the BEADS- 
2D code that is designed to study the dynamics of stellar 
discs in galaxies. The BEADS-2D code is a finite-difference 
numerical code that solves the Boltzmann moment equa- 
tions up to second order in the thin-disc approximation on 
a polar grid (r, 0). More specifically, the BEADS-2D code 
solves for the stellar surface density E, mean radial and tan- 
gential stellar velocities Ur and u^, and stellar velocity dis- 
persion tensor. The latter includes the squared radial and 
tangential stellar velocity dispersions a^^ ^-nd o"^^, respec- 
tively, and the mixed velocity dispersion ct^^. We close the 
system of Boltzmann moment equations by adopting the 
zero-heat-flux approximation. This approximation is equiv- 
alent to assuming the absence of any heat transfer in the 
classic fluid dynamics approach. The BEADS-2D code has 
shown a good agreement with the predictions of linear stabil- 
ity a nalysis of self-gravit ating stellar discs (e.g. Polyachenko 
et al l 19971 . Bertin et aL bOOOD . 

The use of the Boltzmann moment equation approach 
allows us to study non-isotropic stellar systems, in partic- 
ular, the shape and orientation of stellar velocity ellipsoids 
within the discs of spiral galaxies. The BEADS-2D code di- 
rectly evolves observable quantities, and it can follow their 
evolution over many e-folding times starting from arbitrary 
small perturbations. These properties make the stellar hy- 
drodynamics approach advantageous over the more common 
N-body codes which, however, suffer strongly from parti- 
cle noise. Presently, the BEADS-2D code is limited by its 



two-dimensionality (the thin-disc approximation) but a fully 
three-dimensional implementation is currently under devel- 
opment. The interested reader is referred to VT06 for de- 
tails. In the current simulations, the numerical resolution 
has 512 X 512 grid zones that are equally spaced in the (f)- 
direction and logarithmically spaced in the r-direction. The 
inner and outer reflecting boundaries are at 0.2 kpc and 
35 kpc, respectively. We extend the outer computational 
boundary far enough to ensure a good radial resolution in 
the region of interest within the inner 20 kpc. For instance, 
the radial resolution at 1 kpc and 20 kpc is approximately 
10 pc and 150 pc, respectively. 



3 MODEL GALAXY AND INITIAL 
CONDITIONS 

In this paper we study the influence of a spiral gravitational 
fleld on the properties of stellar velocity ellipsoids within the 
disc of a spiral galaxy. Our model galaxy consists of a thin, 
self-gravitating stellar disc embedded in a static dark matter 
halo. The initial surface density of stars is axisymmetric and 
is distributed exponentially according to 



E(r-) = Eo exp(-r/rd). 



(2) 



with a radial scale length of 4 kpc. The central surface 
density Eo is set to 10'^ Mq pc~^. We note that the densities 
of stellar discs are indeed found to decay exponentially with 
distance, with a characteristic scale length increasing from 
2 — 3 kpc for the early type g alaxies to 4 — 5 kpc in the late 
type galaxies (Freeman [l970|). 

The gravitational potential 'I'disc of the stellar disc is 
calculated by 



'I>disc(r, ( 



E(r', <;/>')#' 



r'^ -f r-2 — 2rr' cos( 



(3) 



This sum is calculated using a FFT technique that applies 
the 2D Fourier convolution theorem for polar coordinates 
(Binney & Tremaine IT987I . Sect. 2.8). 

The initial mean rotational (tangential) velocity of stars 
in the disc is chosen according to 



rant 



(4) 



1 



The transition radius between the inner region of rigid ro- 
tation and a flat rotation in the outer part is given by ra^t, 
which we set to 3 kpc. The smoothness of the transition is 
controlled by the parameter Uv, set to 3. The velocity at in- 
flnity, Uao, is set to 208 km s~^. The resulting rotation curve 
is shown in Fig. [l] by the solid line. The initial mean radial 
velocity of the stars Ur is set to zero. 

The radial component of the velocity dispersion is ob- 
tained from the relation arr ~ 3.36(5sGE/«: for a given 
value of the Toomre parameter Qs. Here, k is the epicycle 
frequency. We assume that throughout most of the disc Qs 
is constant and equal 1.3 but is steeply increasing with ra- 
dius at r > 25 kpc. The initial radial distribution of Qs is 
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Figure 1. Rotation curve (solid line) and radial distribution of 
the Toomre Qs parameter (dashed line) in the model stellar disc. 

shown in Fig. [T] by the dashed hne. The tangential compo- 
nent of the velocity dispersion a^^ is determined adopting 
the epicycle approximation, in which the following relat ion 
between cr^^ and a^r holds (Binney & Tremaine Il987l . p. 
125): 

a,, = a..--(^-^ + lj, (5) 

where Uc is the circular speed, i.e. the speed of a hypo- 
thetical star in a circular orbit determined exclusively by 
the gravitational potential. The value of Uc is yet undefined 
(since the gravitational potential of the dark matter halo 
is yet undefined) and we substitute the circular speed with 
the initial rotation velocity of stars u^. We emphasize that 
equation ([5]) is valid only for stars on nearly circular orbits. 
We discuss the practical application of the epicycle approx- 
imation to spiral galaxies in section [5TT] The mixed velocity 
dispersion is initially set to zero. 

Once the rotation curve and the radial profiles of the 
stellar surface density and velocity dispersions are fixed, the 
dark matter halo potential $haio can be derived from the 
following steady-state momentum equation (see eq. (4) in 
VT06) 



dr 



+ 



dr 



-f-— [rT,arr 
r or ^ 



0(6) 



4 DEVELOPMENT OF A SPIRAL 
STRUCTURE 

When starting our simulations we add a small random per- 
turbation to the initially axisymmetric surface density dis- 
tribution of stars. The relative amplitude of the pertur- 
bation in each computational zone is less than or equal 
to 10~^. Since the model stellar disc is characterized by 
Qs ~ 1.3, it becomes vigorously unstable to the growth of 
non-axisymmetric gravitational instabilities. We find that 
the positive density perturbation relative to the initial den- 
sity distribution in the disc can be a good tracer for the 
spiral structure. Hence, we define the relative stellar density 
perturbation as 

E(r,<^)-Si„(r) 



where Ein('') is the initial axisymmetric stellar density dis- 
tribution in the disc and E(r, (/>) is the current stellar density 
distribution. The positive values of ^(r, 0) at f = 1.6 Gyr are 
plotted in Fig.[2]by the red contour lines. The maximum and 
minimum contour levels correspond to the relative pertur- 
bations of 0.8 and 0.1, respectively, and the step is 0.1. The 
position of corotation at approximately 8.5 kpc is shown 
by the dashed circle. The stellar disc has clearly developed 
a two-armed spiral structure and a central bar. A movie 
showing the time development of the spiral structure can 



be found at http:/ /www.univie.ac.at/theis/swing_VT07.avi 



Ei„(r-) 



(7) 



Snapshots from different stages of the evolution have also 
been presented by VT06 for a similar model stellar disc. 

In our numerical simulations, the bar forms within the 
4:1 ultraharmonic resonance (UHR), which is located at ap- 
proximately 4.5 kpc. At the end of the bar, a slightly weaker 
spiral structure emerges almost perpendicularly to the bar. 
Such transitions of bars to spirals are not uncommon (e.g. 
NGC 1672 or NGC 3513). Radially, our spiral arms reach up 
to the OLR 15 kpc). Azimuthally, they cover an angle of 
about TT where they start to fade out. The strongest gradi- 
ents in the stellar density are found near the outer (convex) 
edges of the spiral arms. The pattern speeds of the bar and 
the spiral structure are iden tical. Based on orbital analysis, 
Kaufmann & Contopoulos (|l996l ) identified the UHR as a 
transition region between a bar and spiral features for sys- 
tems characterized by a single patte rn sp eed for both the 
bar and the spiral. Recently, Patsis (|2006l ) has shown that 
such a coupling is related to chaotic orbits that resemble for 
a long time the 4:l-resonance orbital behaviour. We note 
that beyound 10 kpc our spiral arms are very tightly wound 
due to the decreasing pitch angle. This is in contrast to 
logarithmic spirals observed in most spiral galaxies. This in- 
consistency may be caused by the presence of the OLR in 
our model galaxy where spiral arms tend to converge to a 
resonant ring. We also emphasize that the amplitude of the 
spiral stellar density perturbation is quickly decreasing as 
one approaches the OLR, which makes it difficult to observe 
the stellar spirals near the OLR in real galaxies. 

The non-axisymmetric gravitational field of the spiral 
arms and the bar perturbs stellar orbits. To better illustrate 
this effect, we plot in Fig. [5] the residual velocity field of 
stars at t = 1.6 Gyr. The residual velocities are obtained 
by subtracting the initial mean rotation velocities of stars 
(shown in Fig. [T| from the current mean velocities. Notice- 
able large-scale perturbations to the initially axisymmetric 
fiow of stars are seen between the spiral arms and near the 
bar. The fiow pattern near the convex edges of spiral arms 
(and outside corotation) is particularly interesting. The non- 
axisymmetric gravitational field of the spiral arm creates 
two counterpointed large-scale streams of stars, which merge 
near the convex edge of the spiral into a narrow retrograde 
stellar strearrjfl- A similar phenomenon, though less appar- 
ent, is seen near the inner edges of the bar where two merging 



^ We use the notion "stellar stream" to describe a large-scale 
non-circular motion of stars in our numerical simulations. These 
streams may not bo equivalent to the observed moving groups 
of stars, which are characterized by velocities and ages that are 
different from their environment. In principle, multicomponent 
stellar-hydrodynamics or N-body simulations are needed to com- 
pare our stellar streams with the observed moving groups. 
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Figure 2. Contour lines of the positive stellar density perturbation relative to the initial axisymmetric density distribution at t = 1.6 Gyr. 
The maximum and minimum contour levels are 0.8 and 0.1, respectively and the step is 0.1. The residual stellar velocity field in km 
is superimposed on the contour lines (see the text for details). 



streams of stars produce a strongly perturbed velocity pat- 
tern. The regions with merging stellar streams occupy only 
a small portion of the stellar disc. However, as we will see 
in the next section, these regions are characterised by the 
most pronounced influence of the non-axisymmetric gravi- 
tational field on the shape and orientation of stellar velocity 
ellipsoids. 

As demonstrated recently in VT06, the most likely 
physical interpretation for the growth of a spiral structure in 
our model disc is swing amplification. Amplification occurs 
when any leading spiral disturbance unwinds into a trailing 
one due to di fferent ial rotation (Goldreich & Lyndell-Bell 
1 19651 . Toomre Il98ll . Athanas soula Il984 and others). Ac- 
cording to Julian & Toomre (|l966h . the gain of the swing 
amplifier in a stellar disc with Qs — 1.2 and a flat rotation 
curve is the largest when 0.5 <, X <. 2.5. The latter quantity 

is defined as X = A/Acr, where A = 27rr/m is the circum- 
ferential wavelength of an m-armed spiral disturbance and 
Acr = 'in'^G'E/K'^ is the longest unstable wavelength in a 
cold disc. Generally speaking, swing amplification is most 
efficient when A ~ Acr- At A ^ Acr and A <^ Acr, the swing 
amplifier is strongly moderated by local shear and random 
motion of stars, respectively. 

The top panel of Fig. |3] shows the radial profiles of the 
X-parameter for the m — 1 (solid line), m = 2 (dashed line). 



m = 3 (dash-dotted line), and m = 4 (dash-dot-dotted line) 
spiral perturbations. The region of maximum swing ampli- 
fication is bound by the two horizontal dotted lines. It can 
be clearly seen that the m — 1 perturbations are swing am- 
plified in the inner 5 kpc, whereas the m — 2 perturbations 
are swing amplified over a much larger area between 1 kpc 
and 12 kpc. Consequently, the m = 2 perturbations are ex- 
pected to engage more stellar mass into the amplification 
process than the m = 1 perturbations. Hence, it is natu- 
ral to conclude that the m — 2 density perturbations should 
grow faster and ultimately surpass the m = 1 perturbations. 

Why then does our model stellar disc develop a two- 
armed global spiral pattern rather than a three-armed or 
any other multi-armed pattern? Indeed, higher order per- 
turbations are also swing amplified over a considerable area 
of the disc, as Fig. |3] (top panel) demonstrates. The answer 
is that swing amplification needs a feedback mechanism that 
can constantly feed a disc with leading spiral disturbances. 
The trailing short-wavelength disturbances may propagate 
through the centre and emerge on the other side as lead- 
ing ones, thus providing a feedback for the swing amplifier. 
Figure 4 in VT06 and the mentioned movie nicely illustrate 
this phenomenon. However, the propagation of stellar dis- 
turbances through the disc centre is possible if there is no 
inner Lindblad resonance (ILR) or if the waves are reflected 
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before reaching the ILR. A Q-barrier allows for such a re- 
flection of propagating waves (e.g. Athanassoula [l983 ) but 
it is absent in the inner region of our model disc (see Fig.[l|. 

To determine the position of Lindblad resonances 
m(f2p — Q) = ±K, we plot in the bottom panel of Fig. [3] 
the radial profiles of fi and f2 ± w/m for the initial axisym- 
metric stellar disc and two global spiral modes m — 2 and 
m = 3. The angular velocity of the two-armed global spiral 
pattern (see next section) is determined from both visual 
and numerical estimates to be Sip = 23 ± 1 km s~^ kpc~^. 
This value is plotted by the dotted line in the bottom panel 
of Fig. |3] It is clearly seen that the m = 2 disturbances have 
no inner Lindblad resonance, — k/2 < fip at all radii. On 
the other hand, the m = 3 disturbances have the inner Lind- 
blad resonance near the disc centre, f2 — /<;/3 ~ Sip in the 
inner 2 kpc. This prevents the m = 3 (and any higher or- 
der) disturbances from propagating through the disc centre, 
which terminates the feedback mechanism for the swing am- 
plifier for these modes. We find that this mechanism works 
nicely in stellar discs with other radial configurations of the 
X-parameter and pattern angular velocities Qp and can al- 
low for the dominant growth of m = 1, m = 3, and other 
higher order global spiral modes. The results of this study 
will be presented in a future paper. 



5 THE SHAPE AND ORIENTATION OF 
STELLAR VELOCITY ELLIPSOIDS 

5.1 Stellar velocity dispersions and spiral 
structure 

Despite more exotic means of galactic disc heating there 
are currently two known mechanisms involving spirals as 
the heating agents . Thes e are transient spiral stella r density 
waves (e.g Wie len[ l977l . Carlberg fc Se llwood [l985l . Jenkins 

6 Binnev ll990l . De Simone et al. 12004} ) a nd mu ltiple spiral 
stellar density waves (Minchev & Quillen [iooil l . The latter 
mechanism involves two sets of steady state spiral structure 
moving at different pattern speeds. We note that a single 
steady state spiral structure as defined in e.g. Lin et al. 
(|l969h does not heat stellar discs. 

In this section, we focus on the effect that the non- 
axisymmetric gravitational field of a swing amplified spiral 
structure may have on the spatial distribution of stellar ve- 
locity dispersions within a stellar disc. Figures and UJj 
show the spatial distribution of the radial and tangential 
stellar velocity dispersions, respectively, a.t t = 1.6 Gyr. A 
two-armed spiral structure and a bar are apparent in the 
spatial distribution of ffrr and a^^. To better illustrate this 
phenomenon, we plot in Figs |4j; and |4|i the stellar veloc- 
ity dispersion perturbations relative to their initial values 
at t = Gyr. More specifically. Fig. |4j; shows the relative 
perturbations in the radial velocity dispersion calculated as 
(urr— o"rr)/''"rri whcrc a^r is the initial radial velocity disper- 
sion and (Jrr is the model's known radial velocity dispersion 
at t = 1.6 Gyr. To avoid confusions, we note that arr and 
cr^0 are measured in km s~^ throughout the paper. FigureUJi 
shows the relative perturbations in the tangential velocity 
dispersion cr^,/, calculated in a similar manner. The positive 
and negative perturbations are plotted with the shadows of 
red and blue, respectively. The contour lines of the posi- 
tive stellar density perturbations are shown for convenience. 
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Figure 3. Top panel. Radial profiles of the X-parameter for 
the m = 1 (solid line), m = 2 (dashed line), m = 3 (dash-dotted 
line), and m = 4 (dash-dot-dotted line) spiral perturbations. Two 
horizontal dotted lines delimit the region of maximum swing am- 
plification. Bottom panel. Radial profiles of the angular velocity 
of stars n and Q ± K./m for the initial axisymmetric stellar disc 
and two spiral modes m = 2, and m = 3. The angular velocity 
of the spiral pattern Qp is plotted by the dotted line. There is no 
inner Lindblad resonance for the m = 2 mode, whereas the m = 3 
mode has the inner Lindblad resonance near the disc centre. 

Figure |4] clearly demonstrates that the non-axisymmetric 
gravitational field has a pronounced global influence on the 
spatial distribution of stellar velocity dispersions in the disc. 
The most distinctive feature in Fig. |4] is that the positive 
(and negative) relative perturbations in density and both 
stellar velocity dispersions correlate spatially, suggesting a 
causal link. We emphasize that the largest spatial gradients 
in both velocity dispersions (and in the stellar surface den- 
sity) are found near the outer (convex) edges of the spiral 
arms. As we will see in Section [6.11 this may cause a severe 
failure of the epicycle approximation in these regions. 

5.2 The shape of stellar velocity ellipsoids 

In contrast to gaseous discs, stellar discs are essentially non- 
isotropic. This means that the local properties of stars, such 
as mean velocities, are determined by the velocity dispersion 
tensor cr rather than by an isotropic pressure. The tensor cr 
is often non-diagonal in the local coordinate system (r, <j), z). 
The principal axes of a diagonalized velocity dispersion ten- 
sor form an imaginary ellipsoidal surface that is called the 
velocity ellipsoid. The available measurements in the solar 
vicinity indicate that the ratio ai : a2 of the smallest versus 
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Figure 4. Spatial distribution of (a) the radial and (b) tangential 
stellar velocity dispersions at t = 1.6 Gyr. The scale bar is in 
km s^^. The spiral structure is clearly seen in the distribution 
of both velocity dispersions. Perturbations in (c) the radial and 
(d) tangential velocity dispersions at t = 1.6 Gyr relative to their 
initial values. The contour lines of the positive stellar density 
perturbations are shown in (c) and (d) for convenience (see Fig.|2] 
for details). Note that the positive relative perturbations in the 
velocity dispersions correlate spatially with the positive relative 
density perturbations shown in Fig. [2] 



largest principal axes of the stellar velocity ellipsoid in the 
disc plane does not vary signific antly w ith the B — V colour. 
According to Dehnen & Binney (|l998h . the axis ratio cti : (J2 
is approximately 0.6±0.1 for stars with —0.2 ^ B~V 0.7. 
To avoid confusions, we note that throughout the paper cri 
and (72 are measured in km s~^. The B — V colour of an 
ensemble of stars may be considered as a rough estimate of 
their mean age. A weak sensitivity of the axis ratio ai : (J2 to 
the age of stellar populations is intriguing since the principal 
axes of velocity ellipsoids (when considered separately) are 
indeed age sensitive and show a noticeable (a factor of two) 
incre ase along increasing B ~ V colours (Dehnen & Binney 
1 19981 ). 

If the axis ratio ai : G2 is indeed only weakly sensitive 
to the age of stellar populations, it would be interesting to 
determine other phenomena that could possibly affect the 
values of ai : (J2- In this section, we consider the effect that 
the non-axisymmetric gravitational field of the spiral arms 
and the bar may have on the shape of stellar velocity ellip- 
soids within the disc of our model galaxy. Figure [S] shows 
the contour plot of the positive stellar density perturbation 
(red lines) superimposed on the grey-scale map of cri : (72 
in the disc at t = 1.6 Gyr. It is evident that the central 
5 kpc are characterized by ai : (72 ~ 0.8 — 0.9. Most of the 
intermediate and outer disc regions have cri : C72 ~ 0.6 — 0.7, 
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Figure 5. Contour map of the positive stellar density pertur- 
bations (red lines) superimposed on the ratio (Ji : 02 of small- 
est versus largest principal axes of stellar velocity ellipsoids at 
t = 1.6 Gyr. Note abnormally small values of cri : cr2 at the outer 
(convex) edges of spiral arms. 
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Figure 6. Radial profiles of the ratio cri : cr2 (solid line) and pos- 
itive relative perturbation in the stellar surface density (dashed 
line) obtained from Fig.[5]by taking a radial cut at the azimuthal 
angle = 90°. Note that cri : cr2 attains local minima just at the 
outer (convex) edges of spiral arms. 

which is comparabl e to th e solar neighbourhood value of 0.64 
(Dehnen & Binney [19981 ). However, parts of the stellar disc 
that are immediately adjacent to the outer (convex) edges of 
spiral arms are characterized by a noticeably smaller value 
of cTi : (72 « 0.25 — 0.5. To better quantify this phenomenon, 
we calculate the axis ratio ai : o"2 along a radial cut at the 
azimuthal angle (f) = 90° , which is shown in Fig. [5] by the 
dashed line. The resulting ratio cri ; cr2 is plotted in Fig. [S] 
by the solid line, whereas the dashed line shows the relative 
stellar density perturbation C,{r, </>) obtained along the same 
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radial cut. The area below the dashed line is filled with gray 
to better identify the position of spiral arms. The dotted line 
represents the solar neighbourhood value tJi : a2 ~ 0.64. It 
is evident that the axis ratio ai : (72 has deep minima near 
the outer (convex) edges of the spiral arms. The smallest 
value of (71 : (72 ~ 0.25 is attained at r ~ 14.5 kpc. 

The cause for such small values of ai : a-z near the outer 
edges of spiral arms is unclear. One of the principal velocity 
dispersions g\ or (72 becomes much smaller than the other, 
if the mixed velocity dispersion cr^^ becomes comparable to 
the squared radial (p1^^ and tangential ((7^^) velocity disper- 
sions. We will see in the next section that the outer (convex) 
edges of the spiral arms are indeed characterized by large 
values of a^^ and associated vertex deviations. These large 
values of the mixed velocity dispersion are most likely caused 
by the two large-scale (counterpointed) streams of stars that 
merge near the outer edge of the spiral arms into a retro- 
grade narrow stellar stream (see Fig. [5)). It is not clear that 
such streams are a general property of spiral galaxies and 
high-resolution observation are necessary to confirm their 
existence. 

The (possible) existence of regions with abnormally low 
ratios of a\ : G2 can be potentially used to determine the 
position of stellar spiral arms. This task is not trivial since 
the spiral arms in external galaxies are usually manifested 
by strong emission from current star formation, which traces 
the gaseous response to the underlying stellar spiral arms. 
There could exist a substantial phase shift between the 
gaseous and stellar spiral arms. Figure [5] indicates that in 
the inter-arm region (7i : (72 ~ 0.6 — 0.7, while near the 
convex edge of the arm g\ : G2 can become as low as 0.25. 
Such a large difference can be detected. Moreover, our nu- 
merical simulations show that a larger contrast is expected 
for gravitationally more unstable discs. For instance, in a 
Qs = 1.1 stellar disc the g\ : (72 ratio near the convex edges 
of spiral arms can become as small as 0.15 but it stays near 
0.6 — 0.7 in the inter-arm region. On the other hand, multi- 
armed spirals usually show a much smaller contrast in the 
(71 : a2 ratio, possibly due to smaller perturbations to the 
velocity field of stars. High-resolution measurements of the 
stellar velocity ellipsoids in grand-design two-armed spiral 
galaxies are needed to confirm the existence of abnormally 
low values of the axis ratio a\ : (72. 

Finally, it should be stressed that our models are gen- 
uine 2D-models. From a higher-order moment analysis Cud- 
deford & Binney (|l994h showed that neglecting the 2-motion 
can lead to a systematic overestimate of g^^ and, thus, an 
overestimate of the Oort ratio g^^/g^^. Though such 3D- 
effects will definitely have an impact on the exact Oort ra- 
tio, the changes are about 20-25% (cf. e.g. fig. 2 in Cudde- 
ford & Binney). Therefore, we do not expect them to alter 
our results qualitatively, i.e. the large spatial variations of 
(71 : (72 found in our analysis will be kept. In fact, if cr^^ is 
smaller than (7^^, the ratio a'^^/a'^^ should be overestimated 
in our 2D-models and the mentioned variation should be 
even larger. 

5.3 Vertex deviation 

It has long been known from observations that the mixed ve- 
locity dispersion within the disc plane is non-ze ro for stars 
of many stellar types (see e.g. Dehnen & Binnev [19981 ). The 
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Figure 7. Contour map of the positive stellar density perturba- 
tion superimposed onto the map of the mixed velocity dispersion 
(T^^. The scale bar is in (km s^""^)^. 

BEADS-2D code directly evolves the mixed velocity disper- 
sion (7^^ and is particularly suitable for studying the evo- 
lution of (7^^ in galactic discs. Figure [7] shows the map of 
(7^0 at t = 1.6 Gyr. Positive and negative values of g"^^ 
are shown with the shadows of red and of blue, respec- 
tively. We emphasize the fact that (7^^ may be both negative 
and positive. Only the values between —1500 (km s~^)^ and 
-f 1500 (km s~^)^ are plotted to emphasize g^^ in the outer 
disc. It is evident that cr^^ takes large values near the con- 
vex edges of spiral arms and in the inner central region. For 
instance, g^^ can become as large as —1500 (km s"'^)^ near 
the convex edges of spiral arms and ±5000 (km s~^)^ near 
the disc centre. The presence of substantially non-zero C7^^ 
implies that the axis ratio g\ : G2 may become quite differ- 
ent from the ratio g^^ : Grr of the tangential versus radial 
velocity dispersions. 

The non-zero (7^^ indicate that the stellar velocity el- 
lipsoids are not aligned along the local coordinate vector 
r (galactic centre-anticentre direction). The degree of this 
misalignment is often quantified in terms of a vertex devia- 
tion. A class ic defi nition for the vertex deviation is (Binney 
& Merrifield [l998l . p. 630) 

Iv = iatan ( ^ — ^"^^ I • (8) 

According to this definition, the vertex deviation takes the 
values between —45° and +45°. In practice, however, the 
degree of misalignment of the major axis of stellar velocity 
ellipsoids with the centre-anticentre direction may exceed 
these bounds, especially in the regions where the epicy- 
cle approximation breaks down. In an extreme case, one 
might consider a velocity distribution with a vanishing non- 
diagonal element of cr and cr,/,^ exceeding Grr- Then the ma- 
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Figure 8. (a) Temporal evolution of the m = 2 global Fourier 
amplitude A2 (solid line) and mass-weighted vertex deviation Iv 
(dashed line) in our model stellar disc, (b) Dependence of Iv on 
A2 (filled circles). The least-square fit (solid line) yields a near- 
linear relation fv oc A^-'^^ . 



jor principal axis of the velocity ellipsoid will be aligned 
with the azimuthal direction, yielding 1^ = n/2. Therefore, 
we have extended the classic definition (IHll as in VT06 



I. = 



U + sign [cr 



•r A ^ 1 
It cr^r > cr,* 



r<t>} 



— It cr^r 



According to this re-definition, the vertex deviation can take 
the values between —90° and +90°. In this paper we extend 
our previous analysis (VT06) and search for a correlation 
between the magnitude (and sign) of the vertex deviation 
and the characteristics of the spiral gravitational field. 

Global Fourier amplitudes are often used to quantify the 
growth rate of spiral perturbations in gravitationally unsta- 
ble stellar and gaseous discs. These amplitudes are defined 
as 

/Tout /'27r 



1 



E(r, 0, f)e*'"'* r dr i 



(9) 



where m is the spiral mode, Mdisc is the total disc mass, 
and Tin and rout are the disc inner and outer radii, respec- 
tively. The global Fourier amplitudes A{irL) provide a rough 
estimate on the total relative amplitude of a spiral density 
perturbation with mode m. 

The solid line in Fig.|S^ shows the temporal evolution of 
the dominant m — 2 global Fourier amplitude (in log units) 
in our model stellar disc. The m — 2 amplitude indicates 
a steady growth of a two-armed spiral structure in the disc 
from 0.5 Gyr till approximately 1.6 Gyr, when it saturates 
at \ogA2{t) ^ —0.8. The subsequent evolution of the stel- 
lar disc is characterized by a near-constant value of A2{t) 
for approximately another 300 Myr, during which the spi- 
ral structure is being destroyed by the gravitational torques 
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Figure 9. Positive stellar density perturbation superimposed on 
the vertex deviation map at t = 1.6 Gyr. Positive and negative 
vertex deviations are shown with the shadows of red and blue, 
respectively. The numbers indicate the maximum positive vertex 
deviation in the inter-arm region and maximum negative vertex 
deviation in the outer disc. The scale bar is in degrees. 



from spiral arms (this phenomenon is discussed in detail 
in VT06). It is interesting to see the temporal behaviour of 
the vertex deviation as the spiral structure grows, saturates, 
and disperses. The dashed line in Fig. [8^ shows the tempo- 
ral evolution of the mass-weighted vertex deviation (in log 
units) defined as 



E 



dmi /v,i 



(10) 



where dnii is the stellar mass in the i-th computational zone, 
is the vertex deviation in the i-th zone, and summation is 
performed over all computational zones. We plot the mass- 
weighted vertex deviation rather than its arithmetic average 
to deemphasize high vertex deviations in regions of low stel- 
lar density. A correlation between A2{t) and 1^ is clearly 
evident - both grow and saturate in a similar manner. The 
saturated value of the mass- weighted vertex deviation is ap- 
proximately 15°. To emphasize the correlation between A2 
and Iv, we plot in Fig.|8]D the mass- weighted vertex deviation 
versus the m = 2 global Fourier amplitude (filled circles), 
both in log units. The least-square fit shown by the solid 
lines yields a near-linear relation, Iv (x ^42'^^. 

While the mass-weighted vertex deviation appears to 
correlate globally with the amplitude of spiral density per- 
turbations, it is not clear that this correlation should hold 
locally as well. Indeed, Fig. shows the vertex deviations 
in the disc ai t = 1.6 Gyr. The positive and negative ver- 
tex deviations are plotted with the shadows of red and blue, 
respectively. White space corresponds to a near-zero vertex 
deviation. The black contour lines in Fig. |9] delineate the 
spiral pattern by showing the regions with a positive rela- 
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tive stellar density perturbation ^{r,(j)). The minimum and 
maximum contour levels correspond to = 0.1 and = 0.8, 
respectively. It is clearly seen that the value (and sign) of 
the vertex deviation does not correlate with the amplitude 
of the positive stellar density perturbation. In fact, regions 
with maximum amplitude are characterized by near-zero or 
small negative vertex deviations (/v ^ ~5°). At the same 
time, the convex edges of spiral arms, where the stellar den- 
sity perturbations are small but the gradients of the density 
perturbations are large, are characterized by quite large neg- 
ative vertex deviations reaching 1^ = —40° in the outer disc. 
Considerable positive vertex deviations up to -1-30° are seen 
between the spiral arms, where the stellar density pertur- 
bations are negative. Somewhat surprisingly, the bar has a 
near-zero vertex deviation but the central regions on both 
sides of the bar are characterized by the largest vertex devi- 
ations up to /v ~ ±90°. These large values of 1^ are off the 
scale in Fig. [5] (and subsequent figures) to emphasize the 
regions with lower vertex deviation^. 



6 DISCUSSION 

6.1 Comparison with the epicycle approximation 

It is often observationally difficult to measure both radial 
and tangential velocity dispersions of stars. Hence, it is 
tempting to use equation ([Sjl to derive either of the two 
dispersions from the known one and the circular speed of 
stars. One has to keep in mind that equation ([5]) is valid 
only in the epicycle approximation, when the epicycle am- 
plitude (deviation from a purely circular orbit) at a specific 
radial distance r is much less than both r and the scale 
length for changes in both, the stell ar den sity and dispersion 
(Kuijken & Tremaine flOOl . Dehnen [l999l 'l. The latter condi- 
tion may become quite stringent, especially near the convex 
edges of spiral arms. Therefore, it is important to know if 
equation © holds in spiral galaxies with a well-developed 
spiral structure. 

To test the validity of equation ([5]), we calculate the 
ratio : arr from the model's known velocity dispersions 
and compare it with the ratio {a^^ : arr)cp derived using 
equation ([5]). From the observational point of view, it is usu- 
ally difficult to know the exact value of the circular speed 
itc that enters equation ([5]), since it requires independent 
measurements of the axisymmetric gravitational potential 
in a galaxy. The observers usually apply the local and az- 
imuthally averaged rotation velocities of stars and u^, 
respectively, as a proxy for Uc- The validity of this assump- 
tion is checked below. 

To quantify the deviations from the epicycle approxi- 
mation in our model spiral galaxy, we calculate the relative 
errors 

* We note that the sign of the vortex deviation changes when the 
galactic rotation is reversed from the counterclockwise rotation 
(as in present numerical simulations) to the clockwise rotation. 
This happens because the mean tangential velocities (S^ = it^) 
and instantaneous tangential velocities {v^) of stars change the 
sign but the mean and instantaneous radial velocities {vr = Ur 
and Vr, respectively) do not. As a consequence, the mixed velocity 
dispersion defined by equation JTJ and the corresponding vertex 
deviations change the sign, too. 
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Figure 10. Relative errors between the ratio (cr^^ : Urr) de- 
termined from the model's known velocity dispersions and ratio 
((T^^ : (Trr)ep determined from the epicycle approximation ^ 
using the local rotational velocity as a proxy for the circular 
speed itc . The positive and negative relative errors are shown with 
the shadows of blue and red, respectively. The black contour lines 
delineate the positive stellar density perturbation at t = 1.6 Myr. 
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Figure 11. Same as Fig. llUl onlv with the azimuthally averaged 
rotational velocity used as a proxy for the circular speed Uc- 
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using different approximations for the circular speed Uc- In 
particular, Fig. [10] and Fig. Illl show the relative errors that 
are obtained assuming the local rotational velocity and 
azimuthally averaged rotational velocity Urf, as a substitute 
for the circular speed Uc, respectively. The positive and neg- 
ative errors are plotted with the shadows of blue and red, 
respectively. The spiral pattern in the disc is outlined by 
contour lines showing the positive relative stellar density 
perturbations at t = 1.6 Gyr. 

A visual inspection of Fig.[TU]and Fig. Illl indicates that 
the deviations from the epicycle approximation are largest 
near the convex edges of spiral arms. To better illustrate this 
effect, we take a radial cut at (j> — 90° in Fig. 1 101 and Fig. 1111 
(shown by the dashed line) and plot the resulting radial 
distribution of ^{r,<j) = 90°) by the solid lines in Fig. 112b and 
Fig. 112b . respectively. The positive relative stellar density 
perturbation ("(r, (j) = 90°) along the radial cut is shown by 
the dashed lines. The area below these lines is filled with 
grey to identify the position of spiral arms. It is evident 
that if the local rotation velocity of stars is used as a 
proxy for the circular speed Uc fFigure [T2k ). then the relative 
errors can become as large as ^ = +3.0 and ^ = —1.0 near 
the outer (convex) edges of spiral arms. If the azimuthally 
averaged rotational velocity of stars is used as a proxy 
for Uc (FigurefT^Jj), the corresponding relative errors become 
roughly a factor of three smaller. Yet the relative error can 
become as large as ^ = +0.7 near the outer edge of the 
innermost spiral at r ~ 15 kpc, which renders the epicycle 
approximation inapplicable in these regions. 

Finally, we check the validity of the epicycle approxi- 
mation when the circular speed is determined from the ax- 
isymmetric gravitational potential using equation ((6)1 with 
CTrr and ffri,^ Set to zero (cold stellar disc). This approach 
is accurate but difficult to realize observationally for the 
reasons explained above. Nevertheless, it might be instruc- 
tive to compare the ratios of tangential versus radial veloc- 
ity dispersions obtained in our numerical model with those 
obtained from the exact epicycle approximation. The corre- 
sponding relative errors ^(r, are shown in Fig.[l3]with the 
shadows of blue (positive errors) and red (negative errors) . 
The spiral pattern in the disc is outlined by the contour 
lines showing the positive relative stellar density perturba- 
tions at f = 1.6 Gyr. The convex edges of spiral arms are 
again characterized by large relative errors. This effect is il- 
lustrated in Fig. I12fc . which shows ^(r, = 90°) (solid line) 
and C{r,(j) = 90°) (dashed line) calculated along the radial 
cut ai (j) = 90° (Fig. [TS] dashed line). The relative errors can 
become as large as ^ = +0.7 near the convex edge of the in- 
nermost spiral at r 15 kpc and are below ±0.2 elsewhere. 

We conclude that the epicycle approximation is severely 
in error near the convex edges of spiral arms. On the other 
hand, the deviation from the epicycle approximation in the 
inner disc (r < 12 kpc) remains reasonably small, within 
10 — 20 per cent. We argue that one should use the az- 
imuthally averaged mean tangential velocity of stars (if 
available) rather than the local mean tangential velocity u,f, 
as a proxy for Uc and use the epicycle approximation with 
an extreme caution near the spiral arms. It is interesting to 
note that the epicycle approximation holds, to within a 20 
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Figure 12. Radial profiles of the relative error ^{r,<j> = 90°) 
(solid linos) obtained from (a) Figure (TO] (b) Figure [TT] and (c) 
Figure 113! by taking a narrow radial cut at the azimuthal angle 
(p = 90°. The positive relative perturbation in the stellar density 
at t = 1.6 Gyr is shown by the dashed line. 



per cent error, in and near the bar (except for two small 
spots on both sides of the bar where the relative error can 
become as large as ^ = —0.3). 



6.2 Origin of vertex deviation 

It is generally thought that two factors contribute to the ver- 
tex deviation: (i) moving groups and (ii) a large-scale non- 
axisymmetric compo nent o f the galactic gravitational field 
(Binney & Merrifield Il998l l. However, it is not clear which 
of the two factors is more important. The situation may be- 
come even more complicated, since the moving groups may 
partly be caused by the non-axisymmetry in the galactic 
gravitational field. In order to get a better understanding 
of the origin of the vertex deviation in our simulation, we 
analyse the correlation of the vertex deviation with different 
dynamical quantities. 
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Figure 13. Same as Fig. 1101 only with the circular speed deter- 
mined from equation ^ with the velocity dispersions set to zero. 



6.2.1 Correlation of vertex deviation with the 
gravitational field 

We start by examining the role of the non-axisymmetric 
gravitational fiel d in g enerating a vertex deviation. Kui- 
jken & Tremaine ()l994h considered non-axisymmetric grav- 
itational potentials typical for a triaxial (ellipsoidal) dark 
matter halo. They found that the magnitude of the vertex 
deviation should be proportional to the product of sin 2(^b 
and the ratio of the non-axisymmetric versus axisymmetric 
contributions to the total gravitational potential e$. Here 
(j)h is the position angle with respect to the minor axis of 
the equipotential surfaces of an ellipsoid al hal o (see equa- 
tions [14] and [15] in Kuijken & Tremaine [1994 ). The vertex 
deviation is expected to be zero near the minor axis of an 
ellipsoidal halo ~ 0). Of course, a direct comparison of 
our numerical simulations with the analytical predictions of 
Kuijken & Tremaine is problematic, since the dark matter 
halo is axisymmetric in our numerical simulations and the 
non-axisymmetry is introduced by the non-stationary stel- 
lar spiral arms and the bar. Moreover, the simulations al- 
low for a completely non- linear evolution, whereas Kuijken 
& Tremaine's analysis is based on linear perturbation the- 
ory. Nevertheless, it may be helpful to compare two different 
non-axisymmetric structures and their effect on the spatial 
distribution of vertex deviations in a stellar disc. 

In our model, we calculate the non-axisymmetric gravi- 
tational potential by subtracting the gravitational potential 
of the initial axisymmetric stellar disc ^^i^c from the cur- 
rent model's known gravitational potential "l>disc (which is 
the sum of contributions from the axisymmetric and non- 
axisymmetric parts of the stellar disc) . The ratio of the non- 
axisymmetric versus axisymmetric contributions to the total 
gravitational potential in our model is then determined as 
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Figure 14. Contour linos showing the relative perturbations in 
the gravitational potential of the disc superimposed on the spa- 
tial distribution of vertex deviations. The positive and negative 
perturbations are plotted with the solid and dashed lines, respec- 
tively. The positive and negative vertex deviations are shown with 
the shadows of red and blue, respectively. The scale bar is in de- 
grees. 



This quantity calculated at f = 1.6 Gyr is shown in Fig. [14] 
by the contour lines. The minimum and maximum contour 
levels are e# — —0.10 and = -1-0.12, respectively, and 
the step is 0.02. We use dashed and solid lines to differ- 
entiate between the negative and positive values of e#, re- 
spectively. The positive and negative vertex deviations are 
shown with the shadows of red and blue, respectively. It 
is obvious that the magnitude of the vertex deviations in 
our model stellar disc does not correlate spatially with the 
ratio of the non-axisymmetric versus axisymmetric stellar 
gravitational potentials. The maximal vertex deviations (by 
absolute value) in our model are found along the minor axis 
of the central bar, whereas in the case of an elliptical halo 
Kuijken & Tremaine (|l994l ') predicted near-zero vertex devi- 
ations there. The reason for this obvious disagreement is not 
clear and numerical simulations of vertex deviations in mod- 
els with static elliptical dark matter halos will be presented 
in a future paper. 

However, we should keep in mind that the dynamics 
of a star is controlled by the gradient of the gravitational 
potential rather than by the potential itself. We calculate the 
absolute value of the stellar gravitational potential gradient 
in the disc as 

I V$disc| = [id<!>msc/drf + {r~^d<i>disc/d(t>f] . (13) 

This quantity can be regarded as the "strength" of the stellar 
gravitational field in the disc. We further denote the abso- 
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Figure 15. Contour lines showing tlic relative perturbations in 
the absolute value of the gravitational potential gradient of the 
disc superimposed on the spatial distribution of vertex deviations. 
The positive and negative perturbations are plotted with the solid 
and dashed lines, respectively. The positive and negative vertex 
deviations are shown with the shadows of red and blue, respec- 
tively. The scale bar is in degrees. 
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Figure 16. Radial profiles of the vertex deviation by absolute 
value |/v| (solid line) and relative perturbation in the gravitational 
potential gradient x (dashed line) obtained from Fie. llSl bv taking 
a radial cut along the azimuthal angle (f) = 90°. The dotted line 
outlines the position of the bar and spiral arms by showing the 
positive stellar density perturbation. The area below the dotted 
line is filled with grey. 



lute value of the initial stellar gravitational potential gradi- 
ent |V$disc|o and calculate the relative perturbation in the 
stellar gravitational potential gradient as 
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Figure 17. Spatial distribution of (a) derivative dur/dr and (b) 
vertex deviation 1^ in the stellar disc at t=1.6 Gyr. A clear similar- 
ity is seen between the two quantities, especially in the outer parts 
of the disc. Dashed lines show radial cuts at </) = 90° (vertical) 
and ij> = 180° (horizontal). Scale bars are in (a) km s"-'^ kpc"-' 
and (b) in degrees. 



The spatial distribution of </>) at t = 1.6 Gyr is 
shown in Figure [l] with the contour lines. In particular, the 
solid lines delineate positive relative perturbations in the 
stellar gravitational potential gradient and dashed lines de- 
lineate negative relative perturbations. The minimum and 
maximum contour levels are x ~ —0.3 and x = +1-2, re- 
spectively. The positive and negative vertex deviations 1^ 
are shown with the shadows of red and blue, respectively. 
A visual inspection of Fig. [15] reveals a strong correlation 
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between 1^ and positive in the inner disc and a mild 

correlation between these quantities in the outer disc. 

To better illustrate this correlation, we find x along a 
radial cut at (p = 90° (dashed line in Fig. llSp and plot the 
resulting values in Fig. [T5]with the dashed line. The corre- 
sponding vertex deviation is plotted with the solid line. We 
should keep in mind that 1^ changes its sign when the galac- 
tic rotation is reversed but x does not. Hence, we show the 
absolute (rather than actual) value of the vertex deviation. 
The radial position of the spiral arms and the bar at <j) = 90° 
is outlined by the dotted line which represents the positive 
relative perturbation in the stellar surface density (the area 
below this line is filled with grey). The strongest correlation 
between x a-nd |/v| is seen in the central region where both 
quantities attain their maximum values. A noticeable corre- 
lation between |Zv| and x is also seen near the outer (convex) 
edges of spiral arms at r ~ 14 kpc and r ~ 17.5 kpc. On 
the other hand, there is little or no correlation in the inter- 
arm region at r = (15 — 16) kpc, where x is negligible but 
|iv| can become as large as 20°. A similar lack of correlation 
is seen at r = (5 — 10) kpc. We conclude that the relative 
perturbation in the gravitational potential gradient can only 
partly account for the magnitude of the vertex deviation in 
our model disc. 



6.2.2 Correlation of vertex deviation with the streaming 
motions 

In this section, we search for the correlation between the ver- 
tex deviation and local properties of the stellar velocity field, 
in particular, the non-circular streaming motions of stars 
clearly seen in Fig. (2] The streaming motion of stars can be 
quantified by calculating the spatial derivatives of the ra- 
dial and tangential mean stellar velocities, namely, dur/dr, 
dut/i/dcfi, and dur/d(fi. We have examined all three character- 
istics of the streaming motion and found that dur/dr is best 
suited for our purposes. Figure [T7b shows the spatial distri- 
bution of dur/dr (in units of km s~^ kpc~^) in our model 
stellar disc at t = 1.6 Gyr. The positive and negative values 
of this quantity are plotted with the shadows of red and blue, 
respectively. The spatial distribution of the vertex deviation 
at t = 1.6 Gyr, initially shown in Fig. |9l is reproduced in 
Fig. 117b to facilitate the comparison between dur/dr and 
/v A striking similarity between the spatial distributions of 
dur/dr and 1^ is clearly seen in the figure, especially in the 
outer parts of the stellar disc. 

To better illustrate the spatial correlation between 
dUr/dr and l^,, we take two radial cuts at (j> = 90° and 
(j) — 180° (shown with dashed lines in Fig. I17|) and plot 
the resulting radial profiles of dur/dr (dashed lines) and 
/v (solid lines) in Fig. (TH] Now both quantities happen to 
have the same sign and the actual values of the vertex de- 
viation are shown. The comparison of Fig. and Fig. [TS] 
demonstrates that the correlation between dur/dr and 
is in general better than the correlation between x iv. 
Yet there are regions where dur/dr and Iv show no corre- 
lation. For instance, the derivative dur/dr is maximal (by 
negative value) at approximately 14.5 kpc (Fig. 118b .) but 
the corresponding vertex deviation changes its sign and is 
characterized by negligibly small values there. Other char- 
acteristics of the streaming motion {dutf,/d<j) and dur/dcj)) 



show similar correlation with the vertex deviation, though 
less pronounced than in the above case with dur/dr. 

To summarize, the vertex deviation is better correlated 
with the characteristics of stellar streams such as dur/dr 
than with the relative perturbation in the gravitational po- 
tential gradient in the stellar dis(|f|. The relative perturba- 
tion in the gravitational potential shows essentially no cor- 
relation with the vertex deviation, which could be due to the 
fact that the dynamics of stars is controlled by the gradient 
of the gravitational potential rather than the gravitational 
potential itself. It appears that the non-axisymmetric gravi- 
tational field of the spiral arms and the bar perturbs stellar 
orbits by creating large-scale stellar streams. These streams 
introduce (under specific yet poorly understood conditions) 
the mixed velocity dispersion cr^^ into the stellar velocity 
dispersion tensor and consequently generate the vertex de- 
viation. 

This result is in nice qualitative agreement with the lin- 
ear perturbation analysis in KT91. For instance, their equa- 
tion (19) shows also a strong correlation between the vertex 
deviation and the mean velocity field (and its gradients). 
Moreover, the inner bar region is characterized by rigid ro- 
tation, which means KT91's a becomes unity and the de- 
nominator formally vanishes yielding a divergent vertex de- 
viation. Such large vertex deviations are in good agreement 
with the highest values in our simulations found near the 
central bar. 

6.2.3 Higher-order moments 

When applying the Boltzmann moment equations each 
equation of n-th order includes a term of ti 4- 1-order, by 
this formally creating an infinite set of equations. In order 
to end up with a finite set of equations, one has to introduce 
a closure condition. In hydrodynamics this can be a proper 
equation of state (e.g. a polytropic equation of state). In stel- 
lar hydrodynamics the situation is more complex, because 
an obvious equation of state does not exist. For instance, an 
isotropic pressure (or velocity distribution) is guaranteed in 
a "normal" gas by coUisional relaxation which is usually a 
rather fast process compared to other timescales like the dy- 
namical timescale. Contrary to a gaseous system, two-body 
relaxation in a three-dimensional stellar system is a rather 
slow process. Adopting Spitzer & Hart's (|l97lh formula for 
the relaxation time r^ai, 

^3 

" 15.4GV"iln(0.4iV)' ^^^^ 

we get with a velocity dispersion a of 10 kms"'^, a typ- 
ical stellar mass m of 1 Mq, a local mass density p = 
0.2 Mq pc~^ and a number A'^ of starf[f| of 10^, a relaxation 
time of about Trd ~ 1.3 • 10^^ yr. This timescale exceeds the 

^ We stress that the vertex deviation changes its sign when the 
rotation is reversed but the spatial derivatives of the mean stellar 
velocities do not. Hence, for a clockwise rotation, we would have 
obtained an anti-correlation between dur/dr and 1^. 
^ Equation I I15II is derived for a spherical system. Thus we con- 
sider a sphere with a radius equal to a typical disc scale height 
of 100 pc containing about 10^ stars. Though this is only a very 
rough estimate, it does not alter our argument, because N only 
shows up in the Coulomb logarithm term. 
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Figure 18. Radial profiles of the vertex deviation (solid lines) 
and derivative dur/dr (dashed lines) obtained from Fig. 1171 by 
taking radial cuts along the azimuthal angles (a) ij> = 90° and 
(b) cj) = 180° at t = 1.6 Gyr. 

dynamical timescale (and the Hubble time) by far. Even if 
we adopt the smallest values of about 7 kms^'^ for the ve- 
locity dispersions of young stars, the relaxation time is still 
very long. Thus, in a purely stellar galactic disc two-body 
relaxation is negligible, and heat conduction terms which 
are related to the third order moments can be neglected. 

In order to check the third order moments quantita- 
tively we started to perform numerical test particle simu- 
lations in a given time-dependent disc potential including a 
spiral. The numeric al trea tment is similar to the one adopted 
in Blitz fc Sp ergel (|l99lh and more recently in Minchev & 
Quillen ([2003) • First results of these simulations show that 
the 3rd order terms are small compared to the other mo- 
ments, by this confirming our numerical treatment. A more 
detailed analysis of these simulations will be given elsewhere. 



7 SUMMARY AND CONCLUSIONS 

In this paper, we present two-dimensional stellar hydrody- 
namics simulations of a flat galactic disc embedded in a stan- 
dard dark matter halo, yielding a rigid rotation curve in the 
central part and a flat rotation curve outside 3 kpc. We solve 
the Boltzmann moment equations up to second order in the 
thin-disc approximation and follow the growth of a two- 
armed spiral structure deeply into the non-linear regime. 
We analyse in detail the shape and orientation (vertex devi- 
ation) of the stellar velocity ellipsoids when the spiral struc- 
ture has reached a saturation phase. The simulations yield 
the following main results. 

(i) The shape of the stellar velocity ellipsoids (as deflned 



by the ratio ai : 02 of the smallest versus largest axes of the 
stellar velocity ellipsoid) and vertex deviations show large 
variations over the entire stellar disc. In particular, the radial 
distribution of the axis ratio tri : 02 shows deep minima near 
the outer (convex) edges of stellar spiral arms. 

(ii) The spatial distributions of the radial (Jrr and tan- 
gential stellar velocity dispersions show a similar two- 
armed symmetry as the underlying distribution of the stellar 
surface density. 

(iii) The epicycle approximation fails in the regions af- 
fected by stellar spiral arms. In particular, it cannot repro- 
duce the tangential-to-radial stellar velocity dispersion ratio 
04,4, ■ (Jrr obtained from the numerical simulations. The de- 
viations from the epicycle approximation are especially large 
near the outer (convex) edges of stellar spiral arms, where 
the epicycle approximation may be in error by a factor of 
two or more. 

(iv) Globally, there is an excellent correlation between 
the mass-weighted vertex deviation (integrated over the 
whole stellar disc) and the global Fourier amplitude of the 
stellar density perturbations. However, locally, there is no 
simple correlation between the vertex deviation and the 
mass distribution. 

(v) The magnitude of the vertex deviation is not cor- 
related with the gravitational potential (mass distribution) 
and is only weakly correlated with the gradient of the grav- 
itational potential. However, the vertex deviation correlates 
strongly with the spatial gradients of mean stellar velocities, 
in particular, with the radial gradient of the mean radial ve- 
locity, dur/dr 

(vi) Large-scale, non-circular stellar streams induced by 
the non-axisymmetric gravitational field of the stellar spiral 
arms and the bar may be responsible for the failure of the 
epicycle approximation, large vertex deviations, and pecu- 
liar shapes of the stellar velocity ellipsoids. 

Our results are in good qualitative agreement with the 
results of the linear perturbation analysis given in KT91. 
Moreover, our simulations include also the non-linear stages 
of the disc evolution when dynamical feedback limits the 
growth of a spiral structure. The large spatial variations of 
the properties of the stellar velocity distribution found in 
our numerical simulations and the mismatch to the epicycle 
approximation are a caveat for any analysis of observational 
data which cover larger spatial areas or which include stel- 
lar samples coming from a larger spatial region. It will be 
interesting to check whether features like the "u-anomaly" 
could also be created by large-scale streams induced by a 
transient or swing-amplifled spiral. 
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